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We study a lattice model of interacting loops in three dimensions with a l/r^ interaction. Using 
Monte Carlo, we find that the phase diagram contains a line of second-order phase transitions 
between a phase where the loops are gapped and a phase where they proliferate. The correlation 
length exponent and critical conductivity vary continuously along this line. Our model is exactly 
self-dual at a special point on the critical line, which allows us to calculate the critical conductivity 
exactly at this point. 



I. INTRODUCTION 

Continuous phase transitions form a fascinating sub- 
ject in statistical mechanicSfi They are understood es- 
sentially completely in two dimensions^— and can be de- 
scribed in mean field in high enough dimensions. For 
intermediate dimensions, much knowledge is obtained 
from field-theoretic treatments^i^ such as large-iV and 
4 — e expansions, as well as from numerical Monte Carlo 
simulationsi^ In this paper, we consider a class of sta- 
tistical mechanics models with a global U(l) symme- 
try and specific marginally-long-ranged current-current 
interactions2i^ that decay as g/r^ in three dimensions 
(3d), where g is the coupling. Our Monte Carlo study 
suggests that these have continuous phase transitions 
with critical properties such as correlation length expo- 
nent and critical conductivity that vary as a function of 
the coupling g. Our lattice model is self-dual for a spe- 
cial value of the coupling and we know the location of 
the continuous phase transition exactly, a rare instance 
in 3d statistical mechanics problems. This means that we 
know the critical conductivity exactly, but not the cor- 
relation length exponent, for which we need the Monte 
Carlo simulations performed in this paper. 

One of the motivations for our study is the problem 
of phase transitions in matter-gauge systems^^^— with a 
schematic action 



dk in(fc) U,, - ^ ) al{k)a,{k) 



For concretcness, we took bosonic matter fields and also 
wrote the action for the gauge field in fc-spacc. Here we 
will consider singular gauge action with n(fc) ^ \k\/g on 
long wave-lengths in 3d. This would arise, for example, if 
we had a more microscopic model with additional critical 
matter fields^"— that we managed to integrate out (e.g., 
singular gauge field propagators arise in formal large- A'^ 
treatments and the effective coupling g depends on the 
number of flavors). In the present work, we simply pos- 
tulate the singular action at the microscopic level and 



allow ourselves to vary the coupling g at will. We also 
focus on the case of one U(l) matter field. For positive 
m the field ^ is gapped, while for negative m the field 
condenses, and we are interested in the properties of this 
phase transition. 

A convenient representation of the U(l) matter field is 
in terms of integer- valued conserved currents on a lattice, 
with the action 

^matter = ^ ^ " ' ^V) (1) 
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where V^.r. (j' — r') is some short-range interaction. Upon 
integrating out the gauge field, we obtain the action 
in terms of current loops only, with long-range current- 
current interactions. 

A precise definition of a general action for a system of 
loops with long-range interactions is: 



S[J] 



i V{r^r')J,{r)J,{r'). 



(2) 
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Here J are integer-valued currents which reside on the 
links of a cubic lattice with periodic boundary conditions; 
r, r' are the sites of this lattice, and Jfj,{r) is on the link 
between r and r -\- fl. The currents are subject to the 
constraint V- J(r) = for all r. The matter-gauge system 
with the matter action given in Eq. ([1]), upon integrating 
out the gauge field, gives the loop action Eq. ^ with 
the potential V{r - r') = M.r. {r - r') + V^.^. [r - r'), with 
^.r.(fc) = l/n(fc) ~ gl\k\ in fc-space and M.r.(r — r') ^ 
gl\r — r'^ in real space. From now on, we will consider 
such loop-only statistical mechanics models. When the 
overall repulsion in Eq. ^ is large, the system will be in 
a phase with only small loops, while at small repulsion 
the system will be in a phase where the loops proliferate. 



II. MODEL AND MEASUREMENTS 

Before we proceed with the direct Monte Carlo study, 
we first review some results from a duality approach to 



2 



such current loop systems for general V, and will then 
specialize to our precise model. Consider the action 

5orig[Jl ^lY.^{k)\J{k)\^+iY,J*ik)-AUk) . (3) 

k k 

This is the model defined in Eq. ^ written in k- 
space, except that it also includes a static external field 
Aoxt- The partition sum ZfAoxt] can be used to ex- 
tract current correlation functions. We can use a du- 
ality transfor m^^'^^" — to express the partition sum in 
terms of dual variables Q{R), which arc defined on a 
lattice dual to the original lattice. The Q{R) are also 
conserved integer-valued currents, and can be viewed as 
vortex loops ji^i^"— The action becomes 

5dual[Q] = \Y. ^dual(fc) Q{k) + [V X AextKfc) ^ (4) 



Vdual(fc) 



(2^)2 



V{k)\,f{k)\^ 



(5) 



where /^(fc) = f-e*^>, = ^^[2 - 2cos(fc^)] « P 

for small k. This duality is exact for a finite lattice if we 

require Jtot = I]r ^('') = 0' '5tot = 0, which is helpful 
when wc characterize response functions and also when 
we study an exact self-dual point below. 

We monitor loop behavior in our simulations by mea- 
suring the "superfluid stiffness" of the loops, which is 
defined as: 



P''{k) EE 
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where Vol = is the total number of sites. Because 
of the vanishing total current, we measure these at the 
smallest non-zero k. For example, for we used 
fcinin = (0, 0, ^). We can obtain current-current correla- 
tions by differentiating the generating function Z[^ext] 
expressed either in terms of the original variables or 
the dual variables. Equating the two leads to relations 
like^>23 
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for k = (0,0,^2) which is also assumed in the formu- 
lae below. Once again, this relation is exact for a finite 
lattice as long as there are no total current circulations. 

The superfluid stiffness characterizes the current re- 
sponse of the system to an externally applied field A^xt ^ 
However, in a system with long-range interactions there 
is an additional "internal" field created by the system's 
response to Acxt- Thus, if we start with S'mattcr, Eq. ([IJ, 
with fluctuating internal gauge field a governed by the 
action S'gaugo and add the probe field ^oxt, we induce 
non-zero (a). In order to measure the system's response 
to the total fieldj^i^ we define a new observable: 



(k) 
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cr^^(0, Ojfcj) is the Matsubara conductivity^i^^ at the 
imaginary frequency ik^. We can use Eq. ([7]) to show 
that the conductivities in the original and dual variables 
are related by 



'{kK'ik) 



1 



(9) 



In an isotropic system, a^^{k) = cr^^(fc). From now on, 
wc drop Cartesian indices on p and cr. 

We now consider the behavior of p and a in the dif- 
ferent phases of the model. In the small loops phase, 
Eq. ^ gives pjik) ~ /c^, and we can sec from Eq. ([5]) 
that (Jj{k) ^ kz, and so for kz = 2tt/L, both go to 
zero in the thermodynamic limit. In the proliferated 
phase, we know that the dual variables Q are in the 
small loops phase, and we can use the Eq. ([7]) to show 
that I - V{kmin)pjik„un) ^ Kiuai(fcmin)fc2, and Eq. © 

to show Crj(fcmin) ~ l/^z- 

We now specialize to our model, where we take 



V{k) 
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for all k = {kx,ky^kz)] g and t are the parameters of 
the model, with g the strength of the 1/r^ interaction 
and t the strength of an on-site interaction. This model 
was considered by Ref . H, while our study provides a good 
physical characterization and understanding utilizing the 
connection to matter-gauge systems. 

Applying the above analysis, we see that in this 
model l^(fcinin)p(^min) ^ I / L u\ thc Small loops phase 
and ~ 1 — a/L in thc proliferated phase, where a is 
some non-universal number, and we sec that wc can use 
^(fcmin)p(fcmin) m s. manner similar to Binder ratios in 
magnetic systems. Since ^(/cmin) ~ L, we also see that 
p ■ L approaches a finite number in the proliferated phase. 

In our simulations, we will keep g fixed and vary t. 
In addition to p and cr, we measure thermal quantities 
such as average energy (i?), specific heat per site C = 
{{E^) - {Ef)/Yol, and the third cumulant^ 



Vol 
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which can be viewed as a derivative of the specific heat 
with respect to a temperature-like parameter. 

We will also study the derivative of the superfluid stiff- 
ness, dp/dt, which can be evaluated as 



dp{k, t) 
dt 



= (\Uk)\') {E,) - (\Uk)? E,) . (12) 



where Et is the contribution to the energy from thc short- 
range interaction, Et ~ ■ We study thc system 
in the current loop variables J(r) using the directed ge- 
ometric worm algorithm)^^— modified to keep Jtot = 0. 
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icrit(5) 



1/6 0.3125 0.5 2/3 1 1.5 2 



3.2 



3.01 2.50 2.05 1.48 0.98 -1.41 -2.75 -5.69 



TABLE I: Approximate location of our phase transitions, de- 
termined from (7 crossings; tcrit{g) are marked in Fig. [21 Error 
bars are below 0.01. 



III. RESULTS 

At g = 1, the model is exactly self-dual at t = 0, and 
symmetry between J and Q requires that the phase tran- 
sition between the gapped and proliferated phases occurs 
here. For g ^ 1, we must locate the value of t at which 
the phase transition occurs. We did this by holding g 
fixed and sweeping in t. In a model with short-range in- 
teractions, p ■ L is infinite in the proliferated phase and 
zero in the small-loops phase, so one could deduce the 
location of phase transitions from the crossings oi p ■ L 
at different L. We saw above that p ■ L does not diverge 
in the proliferated phase in our long-ranged model, and 
so we are not guaranteed a crossing. However, we argued 
that a is zero in the gapped phase and infinite in the pro- 
liferated phase in the thermodynamic limit, so wc can use 
the crossings of this quantity to determine the location of 
the phase transition. An example of this type of crossing 
is given in Fig. [1] for g = 0.3125. The locations of the 
phase transitions found from these crossings are given in 
Table HI and the phase diagram is shown in Fig. [2] Note 
that our model is ill-defined if any of the V{k) become 
negative. We see from Eq. ((TU)) that this happens first at 
k = (7r,7r,7r), and F(7r,7r, tt) < for t < —gn/^/S. This 
ill-defined region is also labeled in Fig. [51 

At 5 = the model is one of loops with only short- 
ranged interactions, which is well studied, and our iciit 
agrees with existing results <^^i^ For g = 3.2, the sys- 
tem was studied in Ref. |^, and our icrit value agrees with 
theirs. As a further check on our simulations, it is possi- 
ble to derive exact values for (E) and p- L at the self-dual 
point g = 1, t = 0: 



{E){g = l,t = 0;L) 



L3-1 



sin(7r/L) 
2^ ■ 



Our Monte Carlo results are consistent with these rela- 
tions. 

To study the nature of the phase transition, and in 
particular the correlation length exponent we looked at 
the behavior of the derivative L-^. Since p-L approaches 
finite values on both sides of the transition, we know that 
the derivative will be peaked near the critical point. An 
example of the evolution of L • ^ is shown in Fig. [3l 
Finite-size scaling arguments suggest that p ■ L = f[{t — 
tcrit)L^^'^] in our model, and hence the peak value of the 
derivative behaves as 



dt 



(13) 
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FIG. 1: An example of the (T(fcmin) data used to determine 
the location of the phase transition for the model at fixed 
g = 0.3125, varying t. From this data, we determine the 
phase transition to take place at tcrit = 2.05, and estimate 
the crossing value cr"°=='"'5(fcmin) = 0.067 ± 0.005. 
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FIG. 2: The phase diagram of our model, Eq. (|T0|) . The model 
contains a phase with small loops at large t and a proliferated 
phase at small t, separated by a critical line of continuous 
phase transitions. The phase boundary is accurately deter- 
mined from crossings like those in Fig. [Tl (see also Table [l|. 
For large negative t, the model becomes ill-defined. 



The extracted peak values as a function of L are shown 
in Fig. [3^) for a range of parameters. Fitting Eq. 
gives us values of v shown in Fig. [SI 

The same method was applied to C3, which has the 
peak behavior— 



L 



(14) 



with the dimension d — We do not show an example 
of these peaks here, but the extracted [Cajmax values are 
given in Fig. [4)d) and the extracted values of v are shown 
in Fig. [51 Note that although C3 measures a thermal re- 
sponse while p measures a current response, both observ- 
ablcs provide similar values for v. We can see from Fig.[31 
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FIG. 3: L ■ dp/dt as a function of t for the same system as 
Fig. [T] Upon increasing L, the peak evolves slowly towards 
the critical tcrit = 2.05. For each L, we extract the maximum 
value and show it in Fig[4ja). We fit this to Eq. (|13|l to extract 
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FIG. 4: Peak values as a function of L for a) L ■ dp/dt and 
b) Ca, plotted on a log-log scale. We can see that the slopes 
of the lines decrease with increasing g for < <? < 1, which 
corresponds to an increasing v. The Ca data at gr = 2 was not 
shown because it overlaps with the other lines and makes the 
figure hard to read. 



that the curves at large g are not entirely straight on a 
log- log plot. It was unclear what fitting procedure to use 
on this data to extract v, and so we used several different 
procedures. For small g the lines are nearly straight and 
V was not very dependent on the fitting procedure used, 
while for larger g wc observed a larger change in v. The 
error bars in Fig. [5] reflect the different values of v that 
can be obtained by changing the fitting procedure. 

Wc could in principle apply a similar scaling approach 
to the derivative of the conductivity, ^ , but we can see 
from Fig.[T]that these derivatives do not have a maximum 
whose value we can extract. We could instead extract 
values at icrit, but we have found that the results are 
very sensitive to the estimated value of tcrit, and since 
we do not in general know icrit exactly, we were unable 
to obtain precise values of v in this way. At the self-dual 
point, where we do know tent, the extracted v from the 
derivatives at the critical point are consistent with the 
results from the peak values. 

We can see that v changes as we move along the criti- 
cal line in Fig. [21 Though our v values are not deter- 
mined very accurately, they are certainly inconsistent 
with 1/3, which would be the value suggestive of a first- 
order transitions^ In fact, the values for g > are all 
larger than the v of the 3DXY model, and so wc con- 
clude that our phase transitions are second-order. At 
3 = 0, the transition is equivalent to the 3DXY model, 
and our measured value of v is consistent with this. 

For the exactly self-dual model, we were able perform 
simulations at the exact critical point, which allowed us 
to obtain histograms of energy for sizes up to L = 18. 
If the transition were first-order, these histograms would 
have two peaks4^ Our histograms have only one peak, 
and there is no evidence of a 'flat top' which would indi- 
cate a double-peak at larger sizes. This further supports 
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FIG. 5: Values of v determined from scaling of the \L ■ 
cip/cit|max and [Cajmax, (as shown in Fig. 2)) with system size. 



our second-order hypothesis. We were also able to esti- 
mate V using these larger sizes, with results consistent 
with those reported for L < 14. 

During our study of the phase diagram, we obtained 
intersections of curves of a for different sizes. The posi- 
tion of these intersections does not noticeably drift with 
increasing L ioi g < 1, which further suggests a second- 
order transition. In Fig. [6]wc report the values of these 
crossings. The reported values are for L = 12 and 
L ~ 14, since for g > 1 the crossings did drift slightly 
with L. The error bars for these points are a measure of 
this drift, while for points at smaller g the crossings at 
different sizes differ from each other in a non-systematic 
way (presumably due to statistical errors), and the error 
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FIG. 6: The value of cr"°=='"«5 ^g. the couphng g of the 
marginaUy long-ranged interaction. For each g, f,-<=''°=='ng 
determined from plots like those shown in Fig. [T] We can see 
that (j'='^°ss"ig varies as we move along the critical line. The 
crossings are narrow for g < 1, while for p > 1, we noticed a 
weak drift in the crossings as L is increased, and the reported 
values are from the L = 12 and L = 14 data. 



long-range interaction coupling g. We could test this by 
adding additional short-range interactions, e.g. a nearest- 
neighbour interaction, and observing the effect on the 
critical properties. However, we have seen above that 
for the accessible sizes, cr is sensitive to finite-size effects, 
and v has error bars large enough to obscure any small 
change due to the new interaction. We could solve both 
these problems by studying larger sizes. In these simu- 
lations with long-range interactions, increasing the size 
of the system is quite costly at the present time, but 
faster computers or parallelization of the energy calcula- 
tion could make studying these sizes feasible. 

Note that the duality, Eq. inverts the long-range 
coupling g^i^ Due to the short-range interaction, the spe- 
cific model Eq. ([TU]) at g is not exactly dual to the model 
at 1/g. However, if a and ly are dependent only on g and 
not on the short-range coupling, this would imply 



"^(1/5) 



1 



(27r)2o-"°''™s(g)' 



(16) 
(17) 



bars are a measure of these differences. If 0-'^''°^*''"^ is a 
universal quantity, we expect it to be determined by only 
the long-range part of the potential, but for our relatively 
small sizes it is possible that the short-range term makes 
a significant contribution. To get an estimate of this, we 
calculated a without including the short range (t) term 
in the potential in Eq. (jlOp . In the thermodynamic limit, 
this would not change the result. For our data, the mea- 
sured value of was lowered by up to 5%, with a 
greater change observed for smaller g. This error was in 
general smaller than that shown in Fig. [S) 

We conjecture that this intersection point crcrossmg 
function of g only. It varies strongly as we move along the 
line of phase transitions, as shown in Fig. [51 We also see 
that (jcrossmg gyolves continuously into the 3DXY value 
at (7 = 0. From the self-duality of the system at g = 1, 
we know from Eq. Q that 

a-— (fc;.g = l) = ^, (15) 

which is consistent with our data. It is worth noting that 
Eq. (jl5|) is true for all Matsubara frequencies ik^, and in 
particular we can analytically continue to real frequencies 
to obtain the critical dynamical conductivity<^i2& 

IV. DISCUSSION 

We studied a system of loops with 1/r^ interactions, 
and found a line of second-order phase transitions with 
varying critical properties. We were able to exploit the 
duality of the model to get some exact results at 5 = 1. 

It would be interesting to determine whether the crit- 
ical properties of the system are dependent solely on the 



Because v varies slowly and has large error bars, we were 
unable to confirm Eq. ([T7|) . The relation in Eq. ((T5)) is 
satisfied for 0.5 < g < 2, but is not satisfied for the 
pair g = 0.3125 and g — 3.2. However, the crossings at 
5 = 3.2 seem to be drifting towards a value that would 
satisfy the relation as L is increased. Studying the sys- 
tem at larger sizes would enable us to determine whether 
the above relations are satisfied in the thermodynamic 
limit. We could also study the system at larger g, using 
a modified short-ranged interaction to reduce the size of 
the ill-defined region in Fig. [21 

A field-theoretic treatment of the matter-gauge model 
with the singular gauge propagator is possible in the 
spirit of Refs. lol-flll An analysis suggests a line of 
fixed points controlled by the coupling g and contin- 
uously evolving out of the 3DXY fixed point, and it 
would be interesting to study this in detail. It would 
also be interesting to study multi-component systems 
realized with multi-loop lattice modelsi^ii^i^i^"— with 
such marginally-long-ranged interactions. We would 
also like to study models where loops have mutual 
statistics and explore the interplay with the marginal 
interactions jii^^— Loop models with continuously vary- 
ing critical indices can be a fascinating toolbox for study- 
ing phase transitions in three-dimensional statistical me- 
chanics. 
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